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Abstract. The vortex contribution to the probability density function of 
longitudinal magnetization fluctuations is examined in finite 2D XY systems 
close to the Kosterlitz-Thouless-Berezinskii transition temperature. Within the 
temperature range studied their relevance is limited to rare fluctuations, where 
they increase the probability of events exceeding four standard deviations below 
the mean magnetization. 



The characterization of fluctuation statistics is a central problem in the study of 
critical phenomena, as the break down of Landau theory on approaching the critical 
point, implies a non-Gaussian distribution for order parameter fluctuations From 
renormalization group theory it is customary to think of critical phenomena divided 
into universality classes, characterized by the symmetry group of the order parameter 
and the spatial dimension. One would therefore expect critical fluctuation statistics 
to be determined essentially by the set of critical exponents (say (3 and ly for a regular 
critical point) which describe the scaling behavior of derivatives of the singular part 
of the free energy. Evidence from, e.g. the Ising model |j and Potts models Q 
suggests that this is indeed generally the case. 

We have, however recently considered an exception to this established 
phenomenology, the 2D XY model j|, |^. Here the low temperature phase consists 
of a line of critical points, with one independent critical exponent r](T) = T/2i:J 
(with J the exchange constant) extending down to zero temperature and separated 
from the high temperature paramagnetic phase by the Kosterlitz-Thouless-Berezinskii 
(KTB) phase transition [@, |8[ |^. It is well established that for all temperatures 
below the KTB transition temperature, Tktb, renormalisation group flows are to a 
quadratic effective Hamiltonian with the result that the asymptotic behaviour, 

in the limit of large system size is perfectly captured by a harmonic model. The 
advantage of such a simple model is that the probability density function, P(m), for 
fluctuations of the order parameter, m, can be calculated analytically, without using 
either renormalisation group or the scaling hypothesis. Surprisingly, we flnd that 
the form of distribution is independent of along the whole line of critical points. 
This universal scaling function, plotted as the solid line in Fig. 1, arises when aP(m) 
is plotted against /i = (m — (m))/a, where (m) and a are the mean and standard 
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deviation of the distribution jl^, It is asymmetric, with an exponential tail for 

fluctuations below the mean and double exponential for fluctuations above the mean. 

It has been observed that very similar distributions occur in two and three 
dimensional Ising models at a temperature T*{L) slightly below, but close to the 
critical temperature |l^ . We have argued that this is a critical phenomenon [ p^ , 
|l6|, giving weight to our earlier proposition that many correlated systems, both in 
and out of equilibrium, can be driven into a state with very similar fluctuations |l^ , 
irrespectively of their universality class. 

In this Letter we examine the vortex corrections to this universal scaling function 
as one approaches Tktb ^ 0.8929(1) We confirm, on the one hand, that below 

Tktb the effect of vortex pairs will disappears as the thermodynamic limit is taken, 
while above Tktb they become dominant in the same limit. On the other hand, we 
show that their relevance changes only very slowly with length scale, with the result 
that finite size corrections to the thermodynamic limit distribution function will be 
observable over a large range of system sizes. 

The 2D XY model is defined by the Hamiltonian 

= -J^cos(0, (1) 

where the angle 9i gives the orientation of a classical spin vector of unit length, confined 
to a plane (the sum being over nearest neighbors spins) . In the following we consider 
a square lattice of side L with periodic boundaries, and we set throughout J /k-g — 1. 
We define the magnetization m for a single configuration as 

m=l ^ cos(0,-^), (2) 

i=l.N 

where = tan~^ ^i/ Si the instantaneous magnetization direction. 

As the physics of the low temperature phase is perfectly captured by a harmonic, or 
spin wave Hamiltonian |l^, one can, without loss of generality, develop the cosine 
interaction to second order and neglect the periodicity of 9i. This Hamiltonian is 
diagonal in reciprocal space and can be solved straightforwardly. In Refs. [|[ |^ we 
find analytically that aP{m), shown as the solid line in Fig. 1, is a universal scaling 
function, not only of system size, but also of temperature and therefore of critical 
exponent r;. 

We stress that this analytic result is entirely due to harmonic spin waves. At 
any temperature, corrections to it come from two sources: anharnionic, but analytic 
terms in the expansion of Eq. (|l|) and vortex pairs excited as the full periodicity of 9i 
becomes important in the region of Tktb [0 ■ The former have been discussed in detail 
in Ref. |5, 21 ; at fixed temperature, they give small deviations from the universal curve 



for small system size, L. These are finite size corrections to the thermodynamic limit 
function and disappear with relatively modest increase in the system size. 

The two effects can be separated as vortices appear in appreciable numbers in 
a small range of temperature close to Tktb only In order to study the effect 

of vortices on the probability density function (PDF) we have carried out extensive 
Monte Carlo (MC) simulations of the 2D XY model in a small range of temperature 
around Tktb and for different system sizes. We typically use systems of linear size 
L = 16, 32, 64. The system was first equilibrated for 10^ - 10^ MC sweeps (MCs), 
the probability distribution of magnetization was then computed along a trajectory 
of 10^ — 10^ MCs, according to the size of the system. 
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Figure 1. Probability distribution function of magnetization m in a 2D XY 
system of size L = 32, at temperature T around the Tktb transition: T = 
0.86 (□), 0.89 (o), 0.92 (o). In tlie inset the plot is shown in semi-log scale. The 
full line represents the exact solution of the harmonic 2D XY model. 



In Fig. |i| we show results for a fixed system size {L = 32) at several temperatures 
around Tktb- First we observe that, when numerical data are plotted in natural 
units, there is a rather good qualitative agreement with the theoretical curve; a 
more careful inspection however reveals systematic deviations in the tails: when 
observed on a semi-logarithmic scale the exponential tail of the distribution changes 
dramatically, large deviations below the mean magnetization becoming more probable 
below a characteristic breaking point. This is the vortex contribution. The break 
point shifts to higher probability with increasing temperature, consistently with the 
fact that the vortex density increases around Tktb- However at fixed system size 
there is no dramatic difference in behaviour on crossing Tktb. The breaking point is 
a clear signature of two distinct contributions to the PDF ||2^: small fluctuations 
are dominated by spin-waves, while large ones by vortices. It is interesting to 
note, however that, even at Tktb, the vortex dominated region is limited to small 
probabilities and the range |m — (m) | > 4ct, making it difficult to observe. For example 
Zheng and Trimper in Ref. present numerical results concerning an XY system 
of size L = 32 at T = 0.89; but in their plot deviations to the spin wave result do 
not appear as their data is limited to the range \m— {m)\ < Aa. Deviations from 
the spin wave scaling function are also observed in Ref. iQ, however the range of 
temperatures used is outside that for which one expects to see signature of vortices 
for the probabilities resolved. Note that the best agreement with the data originally 
published in Ref. |Q is for the highest temperature studied. 

Let us now consider how does the vortex contribution change as a function of the 
system size. In Figs. 2 and 3 we present the PDF at a fixed temperatures, slightly 
below and above Tktb- In Fig- H with T ~ 0.86 one can see that, although the effect 
of changing size is rather small, the deviation from the spin wave result is reduced 
on increasing L, with the vortex dominated region having a system size dependent 
slope. This is again consistent with a finite size correction to the thermodynamic limit 
function, although it is clear that much bigger sizes would be required to eliminate 
the effects of vortices over the range of probabilities shown here. 
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Figure 2. PDF of magnetization m in a 2D XY system of size L = 
16 (□), 32 (o), 64 (o), at temperature T = 0.86 sligiitly beJow tiie Tktb transition. 
In the inset tiie plot is siiown in semi-log scale (dotted, dash-dotted, and dashed 
lines correspond to L = 16, 32, 64 respectively). 
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Figure 3. PDF of magnetization m in a 2D XY system of size L = 
16 (□), 32 (o), 64 (o), at temperature T = 0.95 slightly above the Tktb transition. 
In the inset the plot is shown in semi-log scale (dotted, dash-dotted, and dashed 
lines correspond to L = 16, 32, 64 respectively). 



^From Fig. || where T = 0.95, one can see that above Tktb just the opposite 
happens; the slope in the vortex dominated region increases slowly with increasing 
L. The region extends out to a cut off, which decreases with system size. The cut 
off appears because, with fluctuations of increasing amplitude, the constraint m > 
comes into play, limiting the range of possible values of fi. For larger system sizes, or 
higher temperature the cut off influences the form of the PDF, the topology changes 
and the system enters the paramagnetic phase through vortex unbinding [ pO| . 

Since the spin-spin correlation length, ^ diverges exponentially with the approach 
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(c « 2), finite size corrections to Tktb are logarithmic in L ||2^, |2^. Defining Tc{L) 
as the temperature where S, = L gives 

Tc(£)^Tk.b+ (4) 
c(logL)^ 

At least within finite size scaling terms, this gives an exceptionally large shift. Our 
results are consistent with this, as at the temperature shown in Fig. ^the correlation 
length is in the range of system sizes studied The results beg the question: can 
one observe scale independence for both the spin wave and the vortex contribution 
to the PDF, giving a universal two component scaling distribution? Data collapse 
onto a single curve along a locus of temperatures with varying system size, would be 
consistent with the notion of the temperature T*{L), introduced in Ref. ||2^ and put 
on a more rigorous footing in Ref. At T*{L), the effective coupling constant, 

renormalized by the vortex pairs in the system, is independent of length scale, over 



all lengths up to size L \2A . It lies in the interval Tktb < T*{L) < Tc{L) and also 
scales logarithmically towards Tktb with system size. It corresponds to a situation 
where the vortex- vortex correlations are independent of scale and one would therefore 
expect two component universality for the PDF at this point. The logarithmic finite 
size scaling means that this is a difficult question to address quantitatively and we 
have not attempted to do so here. However, we do show in Fig. 4 data collected at 
T = 0.93, for three system sizes, much smaller than the thermodynamic correlation 
length. We get fairly good data collapse and dependence on L is extremely small over 
the range studied. There is some deviation from the common curve for L = 64, and 
one could clearly obtain better collapse by moving to a smaller temperature. This 
indicates that at least one locus of points could exist giving data collapse onto a single 
curve. 

In summary, we have presented a quantitative estimation of the vortex 
contribution to the probability density function of longitudinal fluctuations of 
magnetization in the 2D XY model, in a band of temperature above and below 
the Kosterlitz-Thouless-Berezinskii transition temperature. The vortices influence the 
tail of the distribution for large fluctuations below the mean, giving a break in the 
exponential tail towards larger probabilities. In the range of parameter space studied, 
the break is typically for fluctuations of four standard deviations below the mean. 
Given the expected logarithmic flnite size scaling we suggest that this will always be 
the case, within the physical domain probed numerically or experimentally. 
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Figure 4. PDF of magnetization m in a 2D XY system of size L = 
16(n), 32(0), 64 (o), at temperature T = 0.93. Inset: dotted, dash-dotted, and 
dashed hnes correspond to L = 16, 32, 64 respectively. 
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